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£Nj ' At the critical end point in QCD phase diagram, the scalar, vector and entropy suscep- 

t tibilities are known to diverge. The dynamic origin of this divergence is identified within 

the chiral effective models as softening of a hydrodynamic mode of the particle-hole-type 
motion, which is a consequence of the conservation law of the baryon number and the energy. 

o : 

CN ■ §1. Introduction and Summary 

Theoretical studies of QCD thermodynamics have been predicting rich phase 
structure in the plane of the temperature T and the chemical potential 3/i for baryon 
number, as discussed in the workshop, and there are heavy-ion programs now on- 
going and planned in order to explore these new states of matter experimentally. 
Based on the lattice QCD results 1 ) as well as the effective model calculations 2 ) ~~ 5 ) 
it is strongly argued that as one of the prominent features of the real diagram there 
exists a critical end point (CEP)*) of the first-order transition line, which is reviewed 
by M.A. Stephanov 6 ). Since it is a truly singular point, we expect clearer information 
of it from experiments. The basic property of CEP is the divergence of the suscep- 
tibilities of the baryon number and entropy densities as well as of the scalar density. 
Thus it is suggested that the fluctuations in the low-momentum pions and nucleons 
show anomalous behavior as functions of experimental parameters near CEP 7 )' 8 ). 
Theoretical estimates indicate the possibility that the CEP gives still significant ef- 
fects on these observables although the finiteness of the experimental geometry in 
space and time will round off the critical growth of these fluctuations 8 )' 9 ). 

Near the critical point there will be at least one mode which becomes soft repre- 
senting the criticality 10 ) . In this report we address the issue of the soft mode asso- 
ciated with CEP, which causes these critical divergences, within the Nambu-Jona- 
Lasinio model (NJL) 11 )' l2 > and the time-dependent Ginzburg-Landau (TDGL) ap- 
proach 12 ) . The correct identification of the soft mode associated with CEP is crucial 
to determine the time evolution of the critical fluctuations and to study the experi- 
mental signals. 

So far the static properties of CEP have been investigated, where the underly- 
ing (approximate) chiral symmetry is usually emphasized. The effective potential 
with the scalar density a as the order parameter, becomes flat around a nonzero 
equilibrium value of a at CEP, reflecting the divergence of the scalar susceptibility. 
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*' In condensed matter physics the term critical end point is used to indicate a different point 
where a critical line is truncated by meeting a first-order phase boundary. 
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This seemingly suggests the appearance of the gapless sigma mode there. We note 
here, however, that the baryon number and entropy susceptibilities diverge at CEP 
simultaneously. The scalar fluctuation is no more special than others once chiral 
symmetry is explicitly broken*). 

In this report we point out that the most important property of CEP is the 
divergent fluctuations of the conserved quantities such as the baryon number and 
the energy 12 ) . Since only the hydrodynamic mode contributes to the susceptibilities 
of these conserved quantities, the soft mode associated with CEP should have such 
a property and causes various divergences at CEP through the mixing. Within the 
NJL model calculation and the TDGL approach we show that this is indeed the case. 



Divergence of the susceptibility implies that there appears at least one soft mode 
associated with the criticality. Since the susceptibility x 1S obtained as a (/-limit of the 
response function x( u i Q)i which is analytic in the upper half plane of complex to, one 
can express the susceptibility as a sum of the spectral strength over all frequencies: 



where a ultra-violet regularization is understood if necessary. Spectral enhancement 
at vanishing to is required to get the critical divergence because the spectral function 
p(u, q) = 2Imx(w, q) itself is usually integrable. 

Let us discuss the constraints on the spectral functions. First, in the second- 
order chiral transition with the massless quarks, the sigma mode becomes gapless 
so as to realize the manifest chiral symmetry in the spectrum together with the 
pion mode (see Ref. 14) for another possibility). In the typical NJL calculation 15 ), 
the spectral peak of the sigma meson (q = 0) grows sharply and moves to u = 
as the critical point is approached from the symmetric phase, which results in the 
divergent scalar susceptibility. At CEP with the finite quark mass, however, there 
is no symmetry reason to expect the gapless sigma meson. 

The second constraint is a consequence of the conservation of the baryon num- 
ber density and the energy density. The fluctuations of the conserved quantities 
are natural slow modes whose frequencies vanish in the long wavelength limit, and 
constitute a basis of hydrodynamics. For (e.g.,) the quark number density p there is 
a current j such that top — q ■ j = 0, which implies that oj — > as q — > 0. A rephras- 
ing of this fact is that the spectral function of the conserved density fluctuation is 
proportional to the delta function uj5(uj) in the q — > limit 16 )' 12 )' 17 ). We call the 
spectrum with this property "hydrodynamic" in this report. 

At CEP the quark number and entropy susceptibilities are divergent, to which 
only hydrodynamic spectrum contributes. This fact directly requires extra softness 
or critical slowing of a hydrodynamic mode near CEP to give these divergences. 

*' In studying the singularity of CEP, approximate chiral symmetry is at most secondary. The 
size of the critical region of CEP and the distance from TCP are different, interesting issues . 



§2. Susceptibility as a spectral sum and its constraints 
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Fig. 1. Phase diagram of the NJL model in T-/x-m space. Three critical (dashed) lines meet at the 
tricritical point x . The phase boundary (hatched) of the first-order transition forms a wing-like 
structure. 

A remark here may be in order. The coupling with the hydrodynamic mode is 
the origin of the different limiting values of the response function at (u, q) = (0, 0) n ) . 
In fact, the w-limit of x(w, q) 

X (0+ 0) = lim - / Im X (u/, 0) (2-2) 

does not contain the contribution of the hydrodynamic mode, which is apparent from 
a formal relation u8(u) = 0. 

§3. Nambu Jona-Lasinio model n ) 

3.1. Phase diagram 

We study the simple Nambu-Jona-Lasinio model 16 ) C = q(ifl — m)q + g[(qq) 2 + 
{qi^T a q) 2 \ defined with the three-momentum cutoff A and with the coupling strength 
gA 2 = 2.5, within the mean field approximation. The effective potential for the scalar 
density a yields 

f d?k 1 
f2(T, f i,m;o-) = -v j j—^[E k -Tln(l-n + )-Tln(l-n.)] + —a 2 1 (3-1) 

where n± = (e"^") + E k = VM 2 + k 2 , M = m-a and v = 2N f N c = 2-2-3. 

The phase diagram in the T-fx-m space is shown Fig. 1. Within the symmetry 
plane (m = 0), the chiral broken phase is separated from the symmetric phase by a 
critical line, on which the sigma meson mode is gapless. At temperatures below the 
tricritical point (TCP)*) the boundary becomes the first-order line of three-phase 
coexistence 6 ) . As the explicit breaking term m / is introduced and increased, this 
first-order line is lifted to form a wing-like structure. The densities of the scalar, 

Historically the point x is called tricritical because three critical lines meet there. 
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(b) 




Fig. 2. (a) Spectral function in the scalar channel (solid) with \q\/A — 0.1 at a CEP with 
m/A — 0.01. The free gas spectrum (dashed) is also shown for reference, (b) Typical pro- 
cesses contributing to the spectrum. 

baryon number and entropy are discontinuous across this wing surface in general, 
and their susceptibilities contain the critical fluctuations at the edge of the wing*). 

3.2. Spectrum in the scalar channel 

The collective motion of this model is generated through the (pseudo-) scalar 
interactions. The response functions in the scalar and vector channels are easily 
evaluated at the RPA level (a, b = fi,m): 

Xab(iQ4, q) = n ab (iq 4 , q) + n am (iq 4 , q)- 1 .. -2gI7 mb (iq 4 , q), (3-2) 

1 - 2gll mm {iq 4 , q) 

where n ab {iq 4 , q) are the one-loop polarization functions calculated in the imaginary- 
time formalism. The real-time response function is obtained by the usual replace- 
ment iq 4 — > qo + is in the final expression. 

Let us first consider the scalar channel, in which collective mode is generated 
by the denominator in Eq. (3-2). In Fig. 2 (a) we present the spectral function 
p(uj,q) = 21mxm m (uj,q) with \q\/A = 0.1 at a CEP with m/A = 0.01. We clearly 
see a peak around 2M with M/A = 0.331(oc m 1 / 5 ), corresponding to the sigma 
meson, which is interpreted as a collective excitation of the qq pairs (see Fig. 2 
(b)). This sigma meson peak cannot give the critical divergence at CEP because 
it stays massive there. The peak is unphysically narrow in our model because we 
neglected the a <-> hit coupling, which if included will broaden the width and decrease 
the threshold to \J 4m' 2 + q 2 with m n ~ m 2 / 5 / 0. The above conclusion will be 
unaltered in spite of this improvement. 

There is another spectral peak in the spacelike-momentum region {to < |q|), 
which is generated by the collective absorption/emission of the scalar fluctuation 
by a quark or an anti-quark in medium (see Fig. 2 (b)). This type of motion 
is possible only in medium and called particle-hole (p-h) excitation, especially at 
zero temperature. The excitation energy of this motion is kinematically soft; uj ~ 

In the vicinity of the critical lines except for TCP, the mean field approximation is known 
to break down due to large fluctuations. Furthermore, in the chiral broken phase with m = the 
pionic fluctuation is the most important. We neglected these fluctuations in this work. 
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Fig. 3. Spectral function in the spacelike mo- 
mentum region with \q\/A = 0.1, T = T c 
and m/A = 0.01 for several fi (see text). 
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Fig. 4. The pole position of Xmm(<jJ,q) as a 
function of \q\/A for fi/A =0.56 (x), fi c /A 
(*), 0.58 (+). 



-E-fc+q — Eh ~ v ■ q. The peak is strongly enhanced as compared with the free gas 
spectrum shown in the dashed line in Fig. 2 (a). The enhancement of this spectrum 
is demonstrated for n/A=0.58, 0.575, 0.572, 0.571, 0.5701=^/71 (solid lines), 0.5682, 
0.5652, 0.56, 0.55 (dashed lines) with fixed temperature T c /A = 0.19485 in Fig. 3. 

We found numerically a pole of Xmm(^,q) on the negative imaginary axis of cj, 
which we expect is associated with this spectral enhancement. In fact, we confirmed 
that near CEP Xmm(w, q) is well described with this single pole as 

A(g) 

-ijfa+Xmm(q) -iuJ + u c {q) 

with A(q) oc |g| and co c (q) = X(q)x m 1 m (q)- We note that u c {q) oc |g| in the non- 
critical case while u> c (q) oc |q| 3 at CEP because X^mio) ^ Q 2 there (see Fig. 4). 

In the NJL model the singularity of the baryon number susceptibility \nn at CEP 
can be generated again by the denominator in (3-2). It is noteworthy that only the p- 
h collective mode can mix with Xun through the coupling because II m ^(uJ > 0, 0) = 0. 
The importance of the mode with the hydrodynamic character is consistent with the 
argument given in §2. 

We conclude that the origin of the divergences at CEP is the softening of this 
p-h spectrum in the NJL model. We think that the essence of this understanding 
still remain valid after the inclusion of the critical fluctuations beyond the mean field 
approximation although these effects are to be elaborated. If confinement is treated, 
the p-h spectrum may come from the nucleons. 

3.3. Spectral contribution along the critical line 12 ^ 

Now that we know the dominance of the p-h mode in Xmm at CEP, let us 
examine its contribution in the chiral transition with m = because the p-h mode 
always exists in medium. We define ratio R by 

j-, Xmm 

(0,0+) 

Xmm 

R ~ x mm (o,o+) ' (34) 

which measures the importance of the hydrodynamic spectrum as explained in §2. 
This ratio for Xnn is always unity due to the spectral property of the conserved 
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Fig. 5. Ratio R along the critical line. 

quantity. In Fig. 5, we present the ratio R of the scalar channel along the critical 
line approached from the broken phase (see also Fig. 1). It might be a little surprise 
that the p-h mode spectrum gives a finite fraction of the critical divergence even in 
the second-order chiral transition, and it finally dominates the spectral sum (R = 1) 
at TCP, where and xtt diverge in addition to Xmm- This result seems natural 
from the argument given in §2. 

If the critical line is approached from the symmetric phase, on the other hand, 
the p-h collective spectrum does not contribute to Xmm- Thus the divergence solely 
comes from the sigma meson mode (R = 0) and ^ = n /Ifl (0,0 + ). Actually there 
is no p-h spectral strength in the scalar channel at q = in the symmetric phase. 
This is understood as follows: the absorption amplitude of the collective mode with 
momentum q by a quark qi{k) is proportional to a spinor product un{k-\- q)ui{k) oc 
|g| with a helicity flip. (Since a massive quark with a definite chirality has both 
helicity components, the helicity flip is unnecessary in the broken phase.) 

§4. Time dependent Ginzburg Landau approach 12 ) 

In order to confirm the generality of our result, let us repeat our analysis within 
the TDGL approach here. It is essential to introduce another density ip 10 - ) near 
TCP which is physically a linear combination of the baryon number and the energy 
densities, and whose susceptibility diverges at TCP (and at CEP as well). We expand 
the GL effective potential in terms of the ordering densities a and ip around TCP as 

Q = j d 3 x Q(Va) 2 + a a 2 + b a A + ccr 6 + -ya 2 p + ^p 2 -ha- jipj . (4-1) 

Coupling between a and ip must respect the underlying chiral symmetry. Eliminating 
the density p by df2/dip = 7a 2 + p — j = 0, we recover the usual a 6 theory of the 
free energy with a 2 and u 4 coefficients replaced by a = a>o + and b = b — ■ 

In the TDGL approach the coefficients are unknown. In Fig. 6 shown are the 
typical examples of the three-phase coexistence ((a), (b)) and the CEP ((c), (d)). In 
case of exact chiral symmetry (h = 0), the a direction is obviously special. Although 
the GL potential of single order parameter a becomes flat at CEP, the actual flat 
direction is a linear combination of a and p. Both the scalar susceptibility \h = 
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Fig. 6. Schematic GL potentials at a three-phase coexisting point ((a) and (b)) and at a CEP ((c) 
and (d)), with single ordering density (upper) and with two ordering densities (lower). 

l/(2a + 126a 2 +30c<7 ) and the (/9-susceptibility Xj = l + 47 2 <7 2 x/i diverge there. We 
stress here that once h ^ one may equally use an alternative GL potential with 
single order parameter <p defined by eliminating a in favor of ip in Eq. (4-1) 12 \ 

Differently from the NJL case where the collective modes are generated by the 
microscopic interaction and then coupled with the scalar and vector channels, the 
dynamics here in the TDGL approach is introduced phenomenologically from the 
outset 10 ) via 

L a (id t )a = - 6 -^, L v (idt)<p = -^. (4-2) 

One stringent condition on the time dependence is that the density 99 is a conserved 
quantity, and possible form of will be of diffusion L^ui) = —iuj/Xq 2 or sound- 
like Lupiu) = —uj 2 /Xq 2 . On the other hand, it is natural to expect L a (u) = —uj 2 /r 
(propagation) or —itu/r (relaxation). Here we assume coefficients A, r are constant. 

By linearizing the equation of motion around the equilibrium densities, we obtain 
the condition for the normal mode 

d6t ( U{d O^ ^ L ip (id t ) + n w ) =0, (4 ' 3) 

where = 5 2 f2 / 5a5b\ eq . Because of the soft nature of the conserved densities, the 
solution for the mixing between (e.g.,) the propagation and the diffusion is found in 
the small-q 2 limit as 

=f = -(x- h 1 + 4 1 2 a 2 ), j^ = ~Xj 1 . (4-4) 

The soft mode becomes softer in the mixing among the modes, and the normal 
mode with diffusion-like character shows the critical slowing (oc xj 1 —> 0) while the 
propagating mode does not at CEP (a 7^ 0). 
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The response function is the inverse of Eq. (4-3) with the retarded boundary 
condition. Using the spectral function we can again discuss the mode contributions 
to the susceptibilities. For the scalar channel, the response function yields 

Xh[Ul q) " CMw + iF) + + *W - ^ ' 

The susceptibility Xfe = lim ? ^o XkCi is expressed as a spectral sum of the two 
modes as 

Xh = Xh ( ^ X /. 2 2 + .fff 2 2 ) , (4-6) 

where the first term in the bracket is from the propagating mode and the second 
from the diffusion-like mode in Eq. (4-4). From this simple expression we find 
that in the chirally symmetric phase (a = 0) there is no contribution from the 
diffusion mode R = 0. In the second order chiral transition approached from the 
broken phase, a finite fraction of the divergence comes from the diffusion-like mode 
(0 < R < 1), and R = 1 at TCP. In the criticality of CEP the diffusion-like mode 
saturates the divergence (R = 1) because x^ 1 = and a / 0. Incidentally, in the 
u;— limit the diffusion-like contribution drops out due to its hydrodynamic nature. 
These discussions are completely in parallel to the NJL case, and consistent with the 
argument in §2. 
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